Homological Error Correction: Classical and Quantum Codes 
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We prove several theorems characterizing the existence of homological error correction codes both 
classically and quantumly. Not every classical code is homological, but we find a family of classical 
homological codes saturating the Hamming bound. In the quantum case, we show that for non- 
orientable surfaces it is impossible to construct homological codes based on qudits of dimension 
D > 2, while for orientable surfaces with boundaries it is possible to construct them for arbitrary 
dimension D. We give a method to obtain planar homological codes based on the construction of 
quantum codes on compact surfaces without boundaries. We show how the original Shor's 9-qubit 
code can be visualized as a homological quantum code. We study the problem of constructing 
quantum codes with optimal encoding rate. In the particular case of toric codes we construct 
an optimal family and give an explicit proof of its optimality. For homological quantum codes on 
surfaces of arbitrary genus we also construct a family of codes asymptotically attaining the maximum 
possible encoding rate. We provide the tools of homology group theory for graphs embedded on 
surfaces in a self-contained manner. 

PACS numbers: 03.67.-a, 03.67.Lx 
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I. INTRODUCTION 

Quantum Error Correction (QEC) is an important 
breakthrough in the theory of quantum information and 
computation. Without this technique, quantum commu- 
nication over noisy channels would be doom to failure 
and quantum computation would remain in the realm of 
sheer ideal theoretical constructs: powerful in principle, 
but without any chance of being implemented in practice. 

It was Landauer [1], [2], [3] who soon prompted the 
quantum information community to look seriously at the 
problem of quantum errors since they are more harmful 
than classical errors and Unruh pointed out the severe 
negative effects of decoherence [4]. In fact, quantum er- 
rors may show up from different sources: i/ decoherence 
due to undesired coupling of the quantum data with the 
surrounding environment; ii/ imperfections in quantum 
logic gates during the execution of an algorithm. 

The problem of correcting quantum errors seemed 
likely impossible in the beginning, since the classical error 
correcting techniques based on redundancy or repetition 
codes seemed to contradict the quantum no-cloning the- 
orem. Moreover, besides bit-flip errors, there are phase 
errors with no classical counterpart and thus no previous 
theory to compare with. 

Fortunately, all these doubts were dispelled by the 
first quantum error correction code proposed by Shor [5] 
and independently by Steane [6] who showed how to get 
around these difficulties explicitly. Soon, more general 
quantum codes were constructed known as CSS codes [7], 
[8] based on classical correcting codes. These codes are 
very easy to deal with since the correction of bit-flip er- 
rors is factorized out from the correction of phase-flip er- 
rors. CSS codes have found very important applications 
in the security proof of Quantum Cryptography protocols 
without resorting to quantum computers [9]. 

A more general class of codes, encompassing the CSS 
codes, are the stabilizer codes introduced by Gottesman 



[10]. In the stabilizer formalism, the construction of 
quantum codes can be thought of as a task in finite group 
theory for finding Abelian subgroups of the Pauli group, 
leaving invariant a certain subspace which used to en- 
code quantum words. An alternative and independent 
realization was provided by Calderbank et al. [11] using 
the theory of binary vector spaces. 

Despite having a general theory of quantum error cor- 
rection, explicit realization of quantum codes are also 
important in practical implementations. In this regard, 
the number of encoded qubits k, or logical qubits, with 
respect to the number of physical qubits n > k plays an 
important role. The first codes discovered by Shor and 
Steane have a ratio of 1:9 and 1:7, respectively. It is pos- 
sible to show that the best possible ratio for correcting 
one single error is 1:5 [12], [13]. 

The quantum codes mentioned thus far are linear, also 
called additive, codes since the underlying structure is 
that of Abelian stabilizer codes. There are also a series 
of interesting extensions to non-stabilizer codes [14], [15] 
with the aim of increasing the coding capabilities of quan- 
tum codes. For instance, a type of non-additive codes can 
beat the ratio 1:5 of perfect linear codes. It encodes six 
states in five qubits and can correct the erasure of any 
single qubit [16]. A particularly interesting proposal for 
non-abelian quantum codes is due to Ruskai [17] based 
on correcting (2-qubit) Pauli exchange errors besides all 
single qubit errors. This technique can be generalized to 
non- Abelian stabilizer groups based on the permutation 
group S n [18]. 

An alternative approach to quantum error correction 
was introduced by Kitaev [19] known as topological quan- 
tum codes. The notion of topological quantum computa- 
tion was also addressed independently by Freedman [20] . 
This technique allows us to devise topological quantum 
memories which are robust against local errors and ca- 
pable of protecting stored quantum data [21], [22]. 

To understand the notion of a topological code, we 
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first notice that a basic strategy in standard QEC is to 
protect logical qubits by spreading them out in a larger 
set of physical qubits (n > k). This is the reminiscent 
of redundancy in classical codes. In topological quantum 
codes, we go even farther and encode quantum words in 
the nonlocal degrees of freedom of topologically ordered 
physical systems, like certain lattice gauge theories [19], 
[23], [24], [25], or condensed matter systems [26], [27], 
[28], [29], [30], [31]. Detecting topologial order is an im- 
portant issue in this regard [32], [33]. 

Due to this non-local encoding, these quantum code- 
words are intrinsically resistant to the debilitating effects 
of noise, as long as it remains local. This construction 
is rather appealing since it relies on an intrinsic physical 
mechanism for the topological system to self-correcting 
local errors. It means that in a topological code, we do 
not have to check and fix quantum errors from outside 
the system whenever they appear like in standard (non- 
topological) quantum codes. It is the physical properties 
of the system which provide the intrinsic mechanism from 
protecting the encoded quantum states. This mechanism 
is controlled by the interactions described by Hamiltoni- 
ans on certain lattices embedded in surfaces with non- 
trivial topology. The ground state of those Hamiltonians 
exhibit topological order, a type of degeneracy that is ro- 
bust against local perturbations since it is protected by 
a gap from the rest of the spectrum and moreover, the 
degeneracy depends on the topology of the lattice Hamil- 
tonian. Due to this topological order, these states exhibit 
remarkable entanglement properties [34], [35]. 

In addition to being self-correcting, topological quan- 
tum codes exhibit more interesting properties: i/ they 
belong to the class of stabilizer codes; ii/ the interaction 
terms in the Hamiltonian realizing these codes are local, 
i.e., nearest-neighbour interactions. The locality of prop- 
erty ii/ is very important since it facilitates the potential 
physical implementation of these lattice systems. In con- 
trast, the stabilizer operators in non-topological codes 
are generically non-local. 

Practical implementations of topological quantum 
codes have been proposed using optical lattices [36], [37], 
[38] simulating spin interactions in honeycomb lattices 
[26]. In this paper we shall consider only 2-dimensional 
realizations of topological codes, but it is possible to 
make extensions to lattices in 3+1 dimensions [21], [39], 

[40]. ; ; 

The issue of topological quantum computation [19], 
[41], [42], [43], [44], as an instance of fault-tolerance quan- 
tum computing [45], [46], [47], [48], [49], [50], [51] is 
closely related to quantum codes. However, this work 
concentrates only on topological quantum codes. 

In this work we use the terminology of homological 
codes, both classically and quantumly, to highlight the 
fact that they are constructed solely on the information 
about the graph encoded in its homology groups, either 
as simple graphs or as graphs embedded on surfaces. 

The paper is intended to be self-contained and is or- 
ganized as follows: in Sect. II, we introduce the basic 



notions and definitions of classical codes and homology 
groups over Z 2 for graphs. With these tools, we then 
proof theorem II. 2 that allows us to construct classical 
homological codes. Not every classical code is homologi- 
cal, but there exists optimal families of homological codes 
that saturate the classical Hamming bound. In Sect. Ill 
we start recalling the definitions and characterizations 
of quantum codes, then we construct symplectic codes 
for qudits, i.e., quantum states of arbitrary dimension 
D. The idea is to apply the symplectic group ESp£>(n) 
to a trivial code Ct(w, k) of distance one. Symplectic 
codes are equivalent to stabilizer codes. We also intro- 
duce homology of 2-complexes, which are 2-dimcnsional 
generalizations of a graph or 1-complex. With these tools 
we go on to prove theorem III. 2 for constructing qudit 
symplectic codes based on the homology and cohomol- 
ogy groups of graphs embedded in surfaces. Technically, 
these graphs embeddings are called surface 2-complcxes 
that are also introduced earlier. In particular, the cel- 
ebrated Shor's original 9 qubit code can be thought of 
as a homological quantum code belonging to a family of 
codes [{d 2 , 1, d]], with d = 3 (see fig. 20). In general, ho- 
mological quantum codes can be degenerate codes. Next 
we prove a number of important results: 

i/ the subgroup Z 2 appearing in the first homology 
group of non-orientable surfaces is called the torsion sub- 
group. It plays an important role in the construction of 
homological quantum error correcting codes for qudits of 
dimension greater than 2: We show that it is impossible 
to construct these codes with D > 2 on non-orientable 
surfaces, while it is possible to do so for codes based on 
qubits. For orientable surfaces with boundaries, it is pos- 
sible to have homological codes of arbitrary dimension D. 

ii/ we introduce the notion of topological subadditivity 
which is very helpful to find bounds on the efficiency 
(coding rates) of homological quantum codes; 

iii/ for homological quantum codes on the torus, we 
find a family of optimal codes that outperform the orig- 
inal toric codes introduced in [19] and in addition, our 
optimal codes are extended for qudits; 

iv/ we construct an explicit family of quantum homo- 
logical codes for which we can show that the rate k/n 
of logical qubits to physical qubits approaches unity us- 
ing topological graphs embedded on surfaces of arbitrary 
genus; 

v/ it is possible to transform homological codes on 
compact surfaces of arbitrary genus, like the (/-torus, into 
homological codes embedded into planar surfaces with 
boundaries; this is interesting for practical purposes since 
constructing real torus of higher genus does not seem to 
be feasible. 

The results concerning the quantum encoding rate 
were advanced without proof [52] in the particular case 
of qubits (D = 2). 

Sect. IV is devoted to conclusions. In appendix A we 
construct the generators of the symplctic group Spain) 
for the general case of qudits, in appendix B we give a 
detailed explicit proof of the subadditivity property of 
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quantum topological codes, and in appendix C we prove 
that our homological quantum codes for qudits on the 
torus are optimal as far as the coding rate k jn is concern. 

II. HOMOLOGICAL CODES FOR CLASSICAL 
ERROR CORRECTION 

A. Classical error correcting codes 

Classical error correction deals with the problem of 
transmitting messages through noisy channels [53], [54]. 
Usually messages are composed with bits, which can take 
on the values or 1. Such strings of bits, or words, can 
be regarded as vectors over the field Z 2 . The same idea 
holds for the errors introduced in a communication, for 
if u and v are respectively the input and output words, 
we say that the channel has produced the error 

e := v — u. (1) 

An important channel is the (binary) symmetric chan- 
nel. This channel acts on each bit individually, flipping 
its value with certain probability 1 — p, p > ^ . Due to 
the symmetry between and 1, it is possible to assign 
a probability to any given error e, since it does not de- 
pend on the input u. We introduce the weight of a vector 
u £ Z 2 , written wt(u), as the number of non-zero com- 
ponents of u. With this definition, for the symmetric 
channel we have that the probability for a given error e 
to occur is (1 — p) wt ( e ). Thus, errors with small weight 
are more probable, which is important since there is no 
chance to correct an arbitrary error. For words u of in- 
creasing length n, we expect wt(e) ~ np. If we were able 
to correct up to np errors, we would have a successful 
communication with a good probability. 

Given a set of errors S, we say that two words u and 
v are distinguishable with respect to S iff 

Ve,e' e S u + e^v + e'. (2) 

An error correcting code of length n is a subset C of Z 2 . 
Its elements are called codewords. If |C| = 2 fe , we say that 
C encodes k bits. C corrects S if every pair of codewords 
in C is distinguishable with respect to S. Let S(t) consist 
of errors with wt(e) < t. If C corrects S(t) but not S(t + 
1), we say that C is a t-error correcting code. In order to 
characterize this property, let us introduce the distance 
between the words u and v as d(u, v) := wt(u — v). Since 
u+e = v+e' implies d(u—v) = d(e'—e) > d(e')+d(e'), we 
have that two vectors with distance d are distinguishable 
with respect to S(t) iff d > It. The distance of a code is 
the minimum distance between any of its codewords, and 
C is a i-error correcting code iff d > 2t. A code of length 
n, distance d, and encoding k bits is usually denoted by 
[n, k, d]. 

Clearly, the values of n, k and d cannot be arbitrary 
for an [n,k,d] code to exist. In fact, consider a t-error 



correcting code C of length n and |C| = m. Let S n (t) 
contain the elements of S(t) of length n. Since | »SV» (*) j = 
E*=o (?) and u + S n (t)nv + S n (t) = for any pair of 
codewords, we have the (upper) Hamming bound 

m S(") £2 " (3) 

Setting m — 2 k and taking the limit of large n, k, t: 

where n — > as n — > 00 and H(x) is the entropy function 

H(x) := -xlog 2 a;- (1 -x)log 2 (l -x). (5) 

£ is called the rate of the code. A question that nat- 
urally arises here is wether this bound can be reached. 
A theorem by Shannon [55] states that this is asymp- 
totically true, but the codes involved in the proof need 
not be of any practical use. For linear codes, a class 
of codes which we shall introduce below, there is also 
a lower bound known as the Gilbert- Varshamov bound: 
there exists a linear [n, k, d] code provided 




Again, in the limit of large numbers this becomes 

i >(!-*(§)) <!-„), (7) 

where n — > as n — > 00. 

We now focus on linear codes, which have certain prop- 
erties that make them more convenient to use. A lin- 
ear [n, k, d] code is a subspace C of Z^ of dimension 
k for which min uG c-{o} wt(ti) = d. The value for the 
distance follows from the fact that C is closed under 
substraction. A generator matrix G of C is any ma- 
trix with rows giving a basis for C . A parity check 
matrix H for C is any matrix with rows giving a ba- 
sis for C ± , the subspace of vectors orthogonal to any 
vector in C. From this point on, vectors are column 
vectors. To understand why H is useful, first note that 
Hu = u G C. Thus, for any error e and code- 

words u, v we have H(u + e) = H(v + e) = He, that is, H 
measures the error independently of the codeword. He 
is called the error syndrome, and it gives enough infor- 
mation to distinguish among correctable errors. If this 
were not true, then we would have a pair of correctable 
errors such that H(e — e') = => e — e' e C, a contra- 
diction since wt(e - e') < wt(e) + wt(e') < 2t < d. The 
real usefulness of linear codes comes from the fact that 
many codes can be constructed in such a way that the 
deduction of the error from the syndrome is a fast oper- 
ation. As an easy example (due to Hamming), consider 
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FIG. 2: The cycle Ci and the path P 4 . 



FIG. 1: A non-simplicial graph with a self-loop ei and double 
edges e 3 , e 4 . 



the following check matrix for a [7,4,3] code: 



H 



10 10 10 1 
110 11 
1 1 1 1 



(8) 



Notice that columns are the binary representation of 
numbers from one to seven, and thus in this case the 
error syndrome gives the position of the (single) error. 



B. Homology of graphs 

A graph, intuitively, is a collection of vertices and 
edges. Each edge connects two (non necessarily distinct) 
vertices. Figure 1 shows how a graph can be depicted as 
a collection of points or nodes (vertices) linked by curves 
(edges). In such a representation, any intersection of 
edges at points which are not vertices is meaningless. The 
idea of a graph can be formalized in several ways. We 
take here a combinatorial approach, rather than topo- 
logical, and we do not introduce any orientation for the 
edges. 

A (finite) graph T = (V,E,I) (or, if needed, (V r , E T , 
It)) consists of a finite set E of edges, a finite set V of 
vertices and an incidence function / : E — > V(V) such 
that 



l<|/(e)|<2, Vee£. 



(9) 



As usual, Viy) denotes the power set of V, that is, the 
set of subsets of V. The condition over I reflects the fact 
that an edge can only have 1 or 2 endpoints (in the former 
case, it is a self-loop). It is possible to arrange the infor- 
mation conveyed by J in a so-called incidence matrix. To 

\V\ \E\ 

this end, denote V := {vj}- = j and E := {ej}'^. The in- 
cidence matrix has \ V\ rows and \E\ columns. The entry 
in row i and column j is if vi £ I(&j), and 3 — |/(ej)| if 
Vi e I(ej). The incidence matrix for the graph in figure 
1 has incidence matrix 



2 10 
111 
11 



(10) 



Whenever I is not injective we say that T has multiple 
edges. For example, edges e% and e 4 in figure 1 are mul- 
tiple. A graph is called simplicial if it has no self-loops 





FIG. 3: The complete graph K4 and 2 subgraphs, the second 
one a maximal subtree 



nor multiple edges. Note that this is the same as saying 
that the entries of the incidence matrix are or 1 and 
there are no identical columns. 

Two important families of graphs are the n-paths P n 
and the n-cycles C n (n G N). Formally, P n can be 
defined by setting V=l,. . . ,n, E=l,...,n-1 and I(x) = 
{x,x + 1}. For C n , set V = E = Z„ with the same 
description for I. In plain words, n-paths are the com- 
binatorial analog of a closed line segment, while n-cycles 
are the counterpart of a circle. Pictorically, examples are 
shown in figure 2. 

Let 7 and T be graphs. 7 is called a subgraph of T, 
denoted 7 C T, HVy, E~ ( and J 7 are subsets respectively 
of Vr, Er and 7r- We say that two graphs T and T' are 
isomorphic, denoted T ~ T' if there exist two functions 
fj, : Vr — > Vr' and v : Ey — > Ey> which are one-to-one 
and onto and such that 



I r ,(v(e)) = {»(v)\veI T (e)}. 



(11) 



Figure 3 shows some examples of subgraphs. 

A graph isomorphic to some P n is a path, and a graph 
isomorphic to some C n is a cycle. The valence of a ver- 
tex is the sum of the entries in its row in the incidence 
matrix. A path P has one or two distinguished vertices 
with valence distinct of two. We shall call them the end- 
points of P. Two vertices u and v of a graph T are said 
to be connected if there exists a path P C T such that 
the endpoints of P are u and v. This defines an equiv- 
alence relation in V. The equivalence classes are called 
the components of T. A graph with a single component 
is said to be a connected graph. 

A tree is a connected graph with no (sub)cycles. That 
is, a tree is a graph such that for any two vertices there ex- 
ists exactly one path connecting them. Every tree which 
is not a point contains at least two vertex of valence one. 
Some examples of trees are shown in figure 4. A maximal 
subtree of a connected graph T is any tree TCT such 
that Vt = Vr- Equivalcntly, a maximal subtree of T is 
any tree T such that T C T' C T implies that T' is not 
a tree. Thus, there exists a maximal subtree for every 
connected graph. Moreover, given a tree T C T, there 
exists a maximal tree T' such that T C T' C T. 
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FIG. 4: Three trees with one, two and nine vertices. 





FIG. 5: The complete bipartite graph 7^3,3 with a maximal 
subtree (thick line), and the corresponding generating set of 
cycles. 



The Euler characteristic of a graph T, denoted xCO, 
is defined by the formula 



X (r) := |^r| - |£r|. 



(12) 



For any tree T, x(T) = 1 (this can be proved by induc- 
tion on \Vt\)- Thus, if T is any maximal subtree of T, 
then \Er — Et\ = 1 — x(r). For each e G Er — Et we 
define Cr (T, e) as the unique cycle of the graph T + e 
(with the natural definition T + e := (Vr, Et U {e}, It U 
{(e, ir(e))}). The interest of these cycles is that they 
form a maximal set of independent cycles, in a sense that 
will be made clear below. Meanwhile, figure 5 shows an 
example. 

We now introduce the concept of the first homology 
group of a graph T. To this end, we start by defining 
0-chains and 1-chains. Given a graph T = (V,E,I), a 
0-chain is a formal sum of vertices with coefficients in 
Z 2 : 



A„ v, x v e z 2 . 



(13) 



The sum of two chains is defined in a term by term fash- 
ion: 

£A„„ + £A>=£(A„ + A>. (14) 



vev vev 



veV 



We adopt the convention that terms with zero coefficient 
are not written. The especial element with all the coef- 
ficients equal to zero is denoted 0. Let Co(r) be de set 



of 0-chains of T; then (Co(r),+,0) is an abelian group 

isomorphic to zj^' . Note that there is a natural inclusion 
of V in Co giving a basis. The definition of the space of 
1-chains Ci(T) runs along similar lines: just substitute 
V with E. 

Next, we introduce a homomorphism, the boundary 
operator d : C\(T) — > Co(r). It is enough to define its 
value over a set of generators: 



5(e) 




if m 
if m 



{vi,v 2 }, 
{vi}- 



(15) 



It is possible to map naturally subgraphs onto chains; let 
c 7 := X)ee E e i where 7 C Y. Under this identification, 
the boundary of a path with more than one vertex are 
its endpoints, and the boundary of any cycle is 0. 
The first homology group of a graph V is: 



J?i(r) := kerd. 



(16) 



Its elements are always called cycles, but they do not 
necessarily correspond to cycles in the previous sense. To 
avoid confusion, we call the graphs isomorphic to some 
C n simple cycles. We need a description of H\ : 

Proposition II. 1 Let V be a connected graph. Then 
Hi(T) ~ Zl~ x ^ T \ If T is a maximal subtree of V then 
the set { cc T (T,e) \ e & Ep — Et } forms a basis for Hi(T). 
Moreover, if a <E -Hi(r) has coefficients X e on this set of 
edges, then 



Cl = A e C Cr (T.e)- 



(17) 



□ 



If T is composed of several components T t we have: 

^(1^0^). (18) 

i 

Let C°(r) denote the dual space of Cq(T), that is, the 
space of homomorphisms taking Cq(T) into Z 2 : 



C°(r) :=hom(C (r),Z 2 ). 



(19) 



The elements of this space are called 0-cochains. It can be 
regarded as the additive group of functions / : Vr — > Z 2 , 
because a homomorphism is completely defined by giving 
its values on a generating set. Given v E V, we define 

v* e c°(r) by 



v*(u) = S n 



(20) 



where u G V and 5 is the Kronecker symbol. The set 
{v* I v e V r } forms a basis of C°(T). For c° e C°(T), 
c € C (r), we define (c°,c ) := c°(c ). Similarly, C^T) 
denotes the dual space of C\(T) and its elements are 
called 1-cochains. The same comments as for C° are 
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valid substituting V with E, and we use the notation e* 
and (c , Ci) in the same way. 

We define S : C°(T) -> C^r) to be the dual homomor- 
phism of 9, that is, for every c° G C°(T) and ci G Ci(r) 
we have (5c°, c\) :— (c°, dc\). If we think of c° as a func- 
tion over V, then 5c° can be thought of as a derivative or 
gradient. What will be important for us is the fact that 

VveV (5«*,ci) = cieJTi(T). (21) 

If we denote by star(w) the set of edges incident once in 
v, we have 

^* = e *- ( 22 ) 

eEstar(v) 

Although we have maintained our discussion in the 
realm of combinatorics, it is interesting to comment 
briefly how the topological representation of a graph 
T = (V, E, I) is constructed. One starts by giving to V 
the discrete topology. The points of V are called 0-cells. 
We also need a set { D e \ e € E } of closed segments or 
1-cells. The boundary of each of these segments, denoted 
dD ei consists of two points. The information contained 
in / is codified in functions <j) e : dD e — > 1(e) C V with 
the unique requirement that they must be onto. The 
topological space of the graph is then constructed as the 
quotient space of the disjoint union V (J e D e under the 
identifications x <~ 4> e (x) for x G dD e . Properties such as 
connectedness or the first homology group are completely 
topological. 

C. Classical homological codes 

With all the machinery laid down, we are ready to in- 
troduce classical homological error correcting codes. We 
say that a simple cycle isomorphic to C n has length n. 
Let Cy(r) be the set of simple subcycles of T. We in- 
troduce the distance of a graph T, denoted d(T), as the 
minimal length among the elements of Cy(T). 

Given a graph T, let E — {e^}^. Consider the iso- 
morphisms hi : Ci(S) -» Z^ 1 and h 2 : C^S) -» Z^ 1 
defined by 

\E\ 

hi(^A, e4 ) := (Ao,Ai,...,A| B |); (23) 

i=l 

h 2 (^A i e*):=(A ,A 1 ,...,A| B |). (24) 

i=l 

Then 

h 2 (c 1 )-h 1 ( Cl ) = (c 1 ,ci). (25) 

Theorem II. 2 Lef r be a connected simplicial graph, 
not a tree. Construct a parity check matrix H by se- 
lecting a set of linearly independent rows of the incidence 
matrix of T . This gives an [n,k,d\ linear code C with 
n = \E\, k = 1 — x an d d = d(£). 




FIG. 6: The complete graph Ks 

Proof. We claim that hi[i2i(r)] is the code under con- 
sideration. Let F be the subspace generated by the el- 
ements of B := {Sv*\v G V}. From (21) and (25) it 
follows that hi[i?i(r)] = h^F]- 1 . On the other hand, 
since T is simplicial, equation (22) now reads: 

sv* = e *- ( 26 ) 

{eeE\veI(e)} 

Thereby the set of vectors h 2 [B] generates the same space 
as the rows of the parity check matrix H, which proofs 
the claim. 

Since the length is clearly \E\ and k = dim hi [Hi (T)] = 
1 — X, we only have to check the distance of the code. The 
weight function over Z 2 can be pulled back to Ci(r). For 
general 1-chains it gives the number of nonzero coeffi- 
cients in the formal sum. Its restriction to Cy(r) gives 
the length function. Now, let c\ G -ffi(L), c\ ^ 0. There 
exists a subgraph 7 C T such that c 7 = c\ . 7 must con- 
tain a simple subcycle, for if not, then it is a collection 
of trees, and so it contains a vertex v of valence one. But 
then (22) implies (6v*,'y) = 1, a contradiction in view 
of (21). So let c be such a simple subcycle. Clearly, 
wt(hi(ci)) > wt(hi(c)) > d(r), and the equality is ob- 
tained by taking 7 a simple subcycle of minimal length. 

□ 

We do not let T be a tree just to prevent a code encod- 
ing bits of information. Connectedness avoids having 
a code which can be decomposed into two more simple 
ones, but of course there is no problem at all in con- 
sidering unconnected graphs. However, it is completely 
unnecessary to consider a set of disconnected graphs 
Ti since the wedge product of them, \/ i Ti will do the 
work equally well. The wedge product can be obtained 
by choosing one vertex from each graph and identifying 
them all; it does not change the first homology group. Fi- 
nally, if the graph were not simplicial then the distance 
would be 1 or 2, something useless since d — 2t+l, t > 1. 

Let us define u(k, d) as the minimum value of n 
among all the possible [n, k, d] homological codes. Clearly 
v{k,d) < v(k + l,d). In addition, we note that v(k + 
k', d) < v(k, d) + v(k', d), because the wedge product of 
two graphs leading respectively to [n, k,d] and [n',k, d'] 
codes gives a graph associated to a [n + n', k + k', d) code. 
The simplest example of a graph with a code associated 




FIG. 7: The rate k/n vs t/n for the codes generated by the 
families d (□) and K s (A), with the corresponding Ham- 
ming bound (dashed line). Q) is, C-$ = K3. The asymptotic 
Hamming bound is displayed as a reference (continuous line) . 




FIG. 8: The two (connected) topologies for the case k — 2. 
Each curve represents a path, and the labels indicate the num- 
ber of edges composing it. Other topologies are also possible, 
but they could be transformed in one of these by eliminating 
one by one any vertex of valence 1 (an operation which does 
not alter the homology nor the distance). 



is C3. The corresponding code is the repetition code 
{(000), (111)}. This example can be extended to a fam- 
ily of codes in two ways. The easy one is the family Cd of 
[d, l,d] repetition codes. They are clearly optimal, and 
thus, i/(l, d) = d. More interesting is to regard C3 as K 3 . 
In general, the complete graph K s is defined as a simpli- 
cial graph with s vertices and all the possible edges. As 
an example, K 5 is displayed in figure 6. The graph K s 
yields an [(*), (2) — s + 1, 3] code. These codes are clearly 
optimal among homological ones with d = 3. Then we 
can use the family K s to calculate the asymptotical value 
of i/(fc,3). Clearly u(k,3) > k. Let K(s) = (*) - s + 1. 
For k < K{s), v{k, 3) < k) = K{s) + 0(VK(s)). Thus 



lim — = lim — — — - 

t^cot! fc— >oo 3) 



(27) 



in the Ham- 



and asymptotically the point ^ <~ 1, 
ming bound is reached. See figure 7 for a graphical rep- 
resentation of the rates. 

A question that naturally arises is wether every lin- 
ear code is homological. As we shall see, the answer is 
not. Note that the elements of any row of an incidence 
matrix always sum up to two, in Z. So it might be the 
case that a subspace does not have a set of generators 
Ui = (un, . . . ,Ui n ), 1 < i < m fulfilling the condition 



EI 



2 (where the sum must be performed in Z, 



not in Z2). The space generated by the rows of the H 
matrix in (8) is an example of this possibility. To verify 
this, simply check that summing one row to another one 
is equivalent to perform certain column permutation. 

The function f(k, d) behaves well for fixed k = 1 and 
for fixed d = 3. Is this true for other values of the param- 
eters? We do not have a conclusive answer, but a partial 
one may be given. Consider the case k = 2, the (topolog- 
ically) most simple one apart from k = 1. There are only 
two interesting topologies for a graph giving this value 
of k, see figure 8. For case A the inequalities a + b > d, 
a + c > d and b + c > d must hold. Summing up we get 
2n > 3d. The same procedure applied to case B easily 
yields n > 2d. We want n as small as possible, and so in 
principle the first case is the best one. This is confirmed 
by the (optimal) assignment a = b = t + 1, c = t, where 




FIG. 9: A non-homological visualization for the Hamming 
code [7,4,3] . There is one link connecting three nodes, some- 
thing impossible with edges. 



d = 2t + 1. For high values of t, ^ ~ |, and there is 
no way to get a better result. Note how topologies with 
the same first homology group can somehow be classified 
according to their optimality for code composition. If a 
similar calculation is performed for k = 3, K4 is among 
the optimal ones (perhaps as expected) and gives £ ~ \ 
for high values of t. Moreover, due to the high symmetry 
of K s it is possible to construct a bound for its topology 
for any s. One has to consider all the C3 cycles in K s 
and proceed as above to get 



(s-2)n> ( Id. 



(28) 



This is quite a disappointing result, since for high values 
of s one gets £ ~ 0, even for low values of n. However, 
it is not conclusive as long as we do not know wether the 
topology of K s is the optimal one for k — K(s). 

If one does not care about homology and only wants 
a way to visualize codes, then of course it is possible to 
allow "links" connecting an arbitrary number of "nodes" . 
One can still consider chains of links, and the nodes im- 
pose the conditions of the H matrix as vertices did in 
the homological point of view. In such an scheme, the 
code (8) would look as figure 9. It is not clear, however, 
wether this could be of any use. 
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III. HOMOLOGICAL CODES FOR QUANTUM 
ERROR CORRECTION 

A. Quantum error correcting codes 

Quantum error correction is the quantum analogue of 
its classical counterpart. As it usually happens, the quan- 
tum domain gives rise to difficulties not present in the 
classical case; the extension of techniques such as linear 
coding is far from being straightforward. In fact, in the 
early times of quantum information it was believed that 
quantum error correction was impossible. As it happens 
quite often, the dangerous word 'impossible' was soon 
substituted by the more encouraging 'difficult'. 

Here we shall only consider error correction under the 
transmission through quantum noisy channels, which in- 
cludes information storage. This means that we will sup- 
pose that the error correction stage can be performed 
without errors. In general this is not a realistic scenario, 
and the more general framework of fault-tolerant quan- 
tum computation is necessary. 

In the classical case we have been interested exclusively 
in bits. In the quantum case it would seem natural to con- 
sider only qubits, their quantum analogue. Any quantum 
system with only two states can be regarded as a qubit. 
Its representation is a complex Hilbcrt space of dimen- 
sion 2. We shall, however, consider higher dimensional 
quantum systems, or qudits. Although classical compu- 
tation is now far more interested in bits than in dits, it 
was not the case in its early times. In the quantum case 
it seems to be rather interesting to consider qudits [56], 
and thus we will discuss them on equal footing through- 
out this section. 

First a bit of notation. A qudit is described by a 
Hilbert space of dimension D > 2, and finite. This space 
will be denoted V. The elements of a given orthogonal 
basis can be denoted \x) with x — 0, . . . , D — 1. This 
set of numbers is naturally identified with the elements 
of the set modulus 

Z D := Z/DZ. (29) 

In general, whenever an element of Zd appears in an 
expression, any integer in that expression must be un- 
derstood to be mapped to Z^. Messages are strings of 
qudits. Such a string of length n corresponds to the 
space D®". When expressing elements of this space, 
vector notation is useful. As usual, v e ZJ stands for 
v = (vi, . . . , v n ), Vi € Zd- With this notation, we define: 

n 

|v) := (g) k> (30) 

i=l 

The usual scalar product u • v will be employed. It is 
worth noting that, whenever D is not prime, Zjy is not a 
field and Z£, is not a vector space. This is not seriously 
dangerous and we will use the word vector in this wider 
sense [56]. For fixed dimension D, we also introduce the 



symbol 

<p(k):=c^ k (31) 

where k € Z D . 

The essence of quantum error correction is what fol- 
lows. We consider a system S and its environment E. 
The environment cannot be controlled, and it interacts 
with the system producing noise. The system is not ini- 
tially entangled with the environment, but entanglement 
grows with the unavoidable interaction between E and 
S. Omitting the tensor product symbol, this interaction 
can be described as follows [57] : 

|e>|*>->X>fc>M fc |s> (32) 

k 

where |e) and \s) are respectively the initial state of the 
environment and the system, the final states of the envi- 
ronment |e/c) are not necessarily orthogonal or normalized 
and the operators M k acting on the system are unitary. 
In order to perform error correction we need to disen- 
tangle system and environment. This can be achieved 
by enlarging the system S with an ancilla system A and 
whenever it is possible to perform a unitary operation R 
over S' = A ® S such that 

fl(|a)M fc |s» = |a fc >|«> (33) 

where \a) is the initial state of the ancilla. If this is the 
case, then we have: 

J2\ek)R(\a)M k \s)) = fe\e k )\a k )\\8), (34) 
k ^ k ' 

and the errors are gone. Of course, the state \s) is un- 
known. This means that our strategy should work (with 
the same R) for certain subspace of S. Then we could 
use this subspace for information transmission or storage 
without errors. In general, however, we will not be able 
to correct every error and thus we will have to consider 
only errors M k that happen with high probability, just 
as in the classical case. 

Let us explain under which conditions there exists a 
recovery operation R as in (33). A quantum error cor- 
recting code of length n is a subspace C of T>® n such that 
recovery is possible after noise consisting of any combina- 
tion of error operators from some set £ of operators over 
T>® n . The set £ is the set of correctable errors, and we 
say that C corrects £. Note that any linear combination 
of correctable errors is also correctable. A requirement 
for correction to be possible that looks pretty intuitive is 
the following. For every \n) G C such that (£_\rj) = 
and for every M, N e £ 

{€\N*M\rj) = 0. (35) 

This only says that errors do not mix up orthogonal 
states of the code. In what follows we show that in fact 
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this condition is enough and sufficient for (33) to be pos- 
sible. 

Condition (35) can be rewritten in an equivalent form: 
For every j 77) G C and for every N,M G £ 

(Z\NlM\ V ) = c(N^M) (Ch), (36) 

where c(N^M) <G C. Clearly this implies (35). For the 
converse, note that for every |£), \ rj) G C such that (£\rj) = 
condition (35) implies 

0= {t-v\N^M\t + ri) = (Z\NlM\S)-(ri\NlM\Ti), (37) 

from which (36) follows by considering any orthogonal 
basis of C and evaluating N^M on it. 

We now observe that the existence of an ancilla system 
A and a recovery operation R as in (33) implies (36). This 
is because 

{mlMM = ^\M}{^R\a)M 3 \r,) = (^(a^). 

(38) 

The converse is also true; it is enough to take V® n as the 
ancilla system and set 

R(\a)M\0) = M\a)\0 (39) 

for every £ € C and M G £ and for some \a) G C chosen 
as the initial state of the ancilla system. This does not 
define R completely, but it is enough to check that it can 
be extended to a unitary operator over V® n ® V® n . This 
in turn holds true if 

( V \N^M\0(a\a) = (a\N* M\a)( V \$, (40) 

but this follows from (36). 

Our next goal is to introduce a notion of code distance, 
just as in the classical case. A quantum code C is said to 
detect an error N if for every |£), \rj) G C 

(t\N\ V ) = c(N) (fa) (41) 

for some c(N) G C. From the above discussion follows 
that a code C corrects error from £ iff it detects errors 
from the space 

'■= { E NjMi I Mi,Ni e £}. (42) 
1 

For codes of length n, let £ (n, k) be the set of operators 
acting on at most k qudits. We define the distance of the 
code C, denoted d(C), as the smallest number d for which 
the code does not detect £(n,d). Since £(n, ty£(n 7 1) = 
£(n,2t), a code C corrects £(n,t) iff d(C) > 2k. In this 
case we say that C corrects t errors. As in the classical 
case, we can talk about [[n, k, d]] codes when referring to 
codes of length n, dimension D k and distance d. Such a 
code is said to encode k qudits. We use double brackets 
to distinguishing them from classical codes. 

As an example, let us introduce the trivial code of 
length n encoding k qudits: 

C T (n,k):={\0)®\0\(ieV® k } (43) 



where |0) = |0)® n_fe . Since it has distance one, the trivial 
code is quite useless. However, its structure can give rise 
to a rich family of codes. To this end, let U : T)® n — > 
V® n be any unitary operator. Clearly 

UC T (n,k) :={U\c)\ceC T (n,k)} (44) 

is also an error correcting code. In fact, it is clear that for 
any [[n, k, d]] quantum error correcting code C for which 
k is an integer there exists a unitary operator U such that 
U Ct(w, k) = C. These kind of codes are the most usual 
ones. Since 

m\v) = m^UNU^U\r ] ), (45) 

the errors detected by Ct and UCt are in a one to one 
correspondence through conjugation U ■ W. Exploiting 
this idea we could try to find a family of U operators 
for which the calculation of the distance of the code 
U CT{n,k) is easy. This is the subject of the next sec- 
tion. 



B. Symplectic codes 

As a generalization of the usual X and Z Pauli matri- 
ces for qubits, we define for qudits of fixed dimension D 
the operators (31) 

X:= |fc+l><*l. (46) 
kez D 

Z :=J2 ^ k )\ k )^ k V (47) 
kez D 

Note that X D = Z D = 1 and XZ = <p(l)ZX. With 
these operators a basis for the linear operators over V 
can be defined: 

a xz := f(xz)X*Z', (48) 

where x, z G Zd and / : Zn — > C is there to guaran- 
tee a xz = 1. Thus we have to define / by demanding 
f(x) D = tpix) 211 ^ = (-lf^ 1 ), and then we take 

,. . I if D is even and x is odd, . 

/ * == i , -m- aa ■ ■ (49) 

II if L) is odd or x is even. 

The set of c-operators is a basis because 

l*X*l = ^ E rt-W-'Z™. (50) 

ijiGZb 

As an example, note that for qubits we recover the usual 
Pauli matrices: ctoo = I, °"io = A, er i = Z, an = Y. 

We consider strings of qudits. For v G Z 2 ^ and x, z G 
Z£> let us introduce the notation v = (xz) meaning 

v = (xi, ...,x n ,zi,.. .,z n ). (51) 



10 



We can extend our family of operators to act on V® n : But this easily implies that g(v) = w • v for some w e 



»=i 



where v = (xz). We have 



o- x o|v; 

O-0z|v) 



|v + x), 
p(z- v)|v). 



An important commutation relation is [56] 
c u Cv = (p(\i t Qv)a v a u 

where 



1 

-1 



(52) 



(53) 
(54) 



(55) 



(56) 



is a 2n x 2n matrix over Zo- The group of all the opera- 
tors generated by the set of cr-operators is the Pauli group 
Prj(n). Note that there is a natural homomorphism from 
this group onto Z 2 ^ 1 since a u a v oc ct u + v . 

Let us now consider operators U ■ W with U unitary 
such that they are closed over P D (n), that is: 



U(T V U^ = V'(v)cr w ( v ) 



(57) 



Candw: Zg 1 



Z 2 p are functions de- 



where tp : Zj^ 

pending on U. We call this group the extended symplec- 
tic group ESpu(n). It might look that this condition is 
not enough to guarantee that U-W is closed over Pc(n), 
but since it implies ip(v) D = 1, we have ip(v) — ip(g(v)) 



for some g : Z 



In 



Zr>. Thus, there is no problem at 



all. It is can be easily derived that 

C w ( u +v) °C f a (u) ff w(v) °C C7 w(u)+w(v)- 



(58) 



From it this follows that w(u) = Mu where M is a 2n x 
2n matrix over Zp. From (55) we obtain the following 
condition on M: 



M*OM = ft. 



(59) 



The matrix group described by this condition is the sym- 
plcctic group Spa (n) . There is thus a natural group ho- 
momorphism 



h: ESp D (n) -» Sp D (n). 



(60) 



But h is onto, see appendix A, and so it induces the 
isomorphism 



ESpd {n)/ kcr h = Spn (n) ■ 



(61) 



It is interesting to study the kernel of h. For any of its 
elements we have 



Ua v Ui = v{g(v))a v . 



Z 2 r?. On the other hand, 



(63) 



As a result, ker h = Z 2 £ L . 

Now that we have characterized ESpoin), it is time to 
return to our initial purpose of constructing quantum er- 
ror correcting codes. The idea is to apply the symplectic 
group ESpd{ti) to Ct(ti, k) and obtain the codes which 
are called symplectic. A first result is that ker h does not 
help a lot; it only generates codes of the form 



{|u>® |0 ! € e X>® fc } 



(64) 



where u e Z r ^ k . This is an example of conjugated codes. 
More generally, for each symplectic [[n, k, d\] code C there 
exists a family of D n ~ k conjugated [[n, k, d\] codes ob- 
tained from C by application of cr-operators. This will 
become clear shortly. As a result, we only have to focus 
on Spi)(n) when looking for better codes. 

For any subspace V C Zj^ we define the subspace 



V := { u e Z 2 S | Vv e V u*fiv = }. 



(65) 



If V C V we say that V is isotropic. Now let Vq t C Z 2 ^ 
be the isotropic subspace containing the elements of the 
form (Oz), where z € Z^, must have its last k elements 
equal to zero. It is not difficult to verify that Cxin.k) 
detects tr v iff 



v £ V Ct -V Ct . 



(66) 



Consider any symplectic code C = U Crin^k) with 
h(U) = M. Wc can define V c := MV Ct , giving V c = 
MV Ct = MV Ct . Then C detects a v iff v g % - V e . In 
analogy with the weight function for classical codes, for 
any v = (xz) G Zg\ let 

|v| := \{i = l,...,n\ Xi ^0 or m ^0}| (67) 

Recall that cr-operators over one qudit form a basis. This, 
the fact that the space of operators detected by a code 
is a linear subspace and the previous discussion imply 
altogether: 



d(C) = min |v|. 

vEVc-Vc 



(68) 



This equation shows that the distance of the code de- 
pends only upon Vc- On the other hand, given two 
isotropic subspaces V\,Vi C Z 2 £ of the same dimension 
it is possible to find a matrix M <G Spo(n) such that 
MV\ — V2 [56]. Therefore, for any isotropic subspace of 
dimension n — k such that V C V there exists an [[n, k, d]] 
symplectic code C with Vc = V. This way, the problem 
of finding good codes is reduced to the problem of finding 
good isotropic subspaces V C Z|j\ This is analogous to 



(62) the classical situation with linear codes. 
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It is worth revisiting the trivial code on a new light. 
Consider the following abclian subgroup of P/j (n) : 



S T (n, k) :— { cr v | v e V Ct }■ 



(69) 



The trivial code can be defined just in terms of this group: 

C r (n, k) = { |0 e P 8 " |Vcre 5r(n, k) a\£) = |£) }. 

(70) 

Sr{n,k) is called the stabilizer of Cr(n, fc). The stabi- 
lizer of any code C = U CT(n,k) is the abelian group 
Sc ■= U ST(n 7 k)U^ , and C can be defined by its sta- 
bilizer just as we did for Cr(n, fc). It is because of this 
point of view that symplcctic codes are also called stabi- 
lizer codes. A question that naturally arises here is under 
which conditions an abelian subgroup S C Pr>(n) is the 
stabilizer of a symplcctic code. Clearly S must fulfill the 
condition 



(71) 



For D prime this is the end of the story, but in other 
case a bit of care is necessary, as we shall show now. 
Because of condition (71), S is isomorphically mapped to 
a subgroup Vs C Z 2 £ . We claim that S is the stabilizer 
of a symplcctic code iff Vs is a subspace of Zj^ [56] . We 
only have to check the if direction. First, the elements 
of S can be labelled with the elements of Vs- We denote 
them o\s(v), v € V5. Vs is isotropic, and so we can find a 
symplectic code C such that Vc — Vs- Let us denote the 
elements of its stabilizer ctc(v), but in such a way that 



0c(v) = <p{g(v))a s {v), 



(72) 



where g : Zf} — ► C and v e Vc = V5. This is always 
possible since <Js(v) D = 0"c( v ) D = 1- It is easily verified 
that g is linear, but then g(v) = w • v for some w € Z^ 1 . 
Due to (63), there is a conjugate code of C such that S 
is its stabilizer. 

Although condition (35) guarantees that recovery is 
possible, it is worth giving a more concrete recipe for 
symplcctic codes. So let C be a code of distance d, {vi} 
a basis of Vc and G := {<f{fi)<J Vi } a generating set for 
its stabilizer, where /j <E Zpj- Suppose that an encoded 
state |£) has been subject to correctable noise as in (32): 



k 



|ek>0- u J£}. 



(73) 



where <G Zj^ and |ufc| < d/2. We first measure the 
syndrome of the error. This amounts to project the sys- 
tem to any of the eigenstates of each operator in G. For 
each of the eigenstates there is a corresponding eigenvalue 
ip(gi), g% <E Zd- The final state then is proportional to 



|e*>u»IO- 



(74) 



Let <7 U and ov be any of the error operators in this sum. 
Note that u— u' e Vc- Also, ctLct u oc o- u -u' is detectable, 



and so in fact u — u' e Vc- With the information from the 
error syndrome, we can choose any w such that u*f2w = 
gi and ct w is correctable. Then any of the error operators 
in the sum is of the form o~ Uk oc <7 W (J U ' with u' k £ Vc- In 
other words, (74) can be rewritten 



14) 0wia 



(75) 



where \e' k ) oc |e^). This means that the measurement by 
itself is enough to disentangle system and environment, 
and we only have to perform to recover the original 
encoded state. 

Due to the essential role of Vc, symplectic codes are 
usually given in the form of a 2n x (n — k) matrix whose 
rows form a basis for it. As an example, there is a 
[[5, 1,3]] symplectic code [12], [13] of the form 



1 1 1 1 

1 1 1 1 

1 1 1 1 

1 1 1 1 



(76) 



An important class of codes is that of the so-called CSS 
codes. For this codes the matrix has the form 



H 
H 



where H is the check matrix of a classical code C such 
that C cC 1 . In fact, more generally, any code for which 
the matrix can be put in the form 



H x 
H 2 



in such a way that X and Z operators are not mixed up, 
is called CSS. 

Returning to general codes, it is possible to derive a 
quantum analogue of the Hamming bound for certain 
quantum codes. Let 



£ a {n,t) := {<7 V I |v| <t}. 



(77) 



It is clear that a code that corrects £ a (n,t) corrects t 
errors. Let C be a code of length n and dimension m 
that corrects t errors and satisfies the condition that for 
every normalized |£) e C and for every u, v e Z 2 ^ such 
that |u|, |v| < t 



(78) 



Such codes are called orthogonal or nondegenerate. No- 
tice that for V® n there are (D 2 - 1)*(") er-opcrators of 
weight t. This and condition (78) give the quantum Ham- 
ming bound [58] 



E(^ 2 



iy\ ,v ) < D n . 



(79) 
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FIG. 10: A 2-complex composed of 9 vertices, 21 edges and 4 
faces. 



C. Homology of 2-complexes 

A 2-complex is the 2-dimcnsional generalization of a 
graph or 1-complex. In general one can speak of cell 
complexes of arbitrary dimension, but we will keep things 
simple and restrict our attention to these low-dimensional 
cases. Recall that graphs were obtained by attaching 
1-cells (arcs) to a set of 0-cells (points). We can con- 
tinue the process by attaching 2-cells (discs) to the graph. 
Here attaching means "identify points in the boundary 
through continuous maps"; recall the end of section II B. 
Indeed, we will not consider such general 2-complexes. 
We are interested in the combinatorial point of view, and 
our definition will reflect this fact. Figure 10 shows an 
example of the kind of objects we shall consider. The 
goal is to study the first homology group Hi of these 
objects. Although our study of graphs only included Z2 
homology, now we will discuss Z homology. In fact, when 
we talk about qudits we will be interested in Zd homol- 
ogy, but this is constructed substituting Z for Z D in the 
definitions. 

Moving from Z2 homology to Z homology requires the 
introduction of orientation. An oriented finite graph T — 
(V, E, I s , It) consists of a finite set V of vertices, a finite 
set E of edges and two incidence functions I s , I t : E — ► 
V. The subindexes stand for 'source' and 'target'. We 
say that an edge e € E goes or points from I s (e) to 
It(e). Let us introduce the set of inverse edges E^ 1 := 
{e _1 I e G E}, where e _1 is just a symbol and we set 
(e^ 1 )- 1 := e. We will use the notation E := E U E' 1 . 
The incidence functions can be extended to E setting 
lg(e) = Itie" 1 ) for any e e E. 

In order to give a combinatorial meaning to the attach- 
ment of discs to graphs described above, we introduce the 
idea of walks on graphs. Given an n-tuple (a, 6, c, . . . ), 
let [a, b, c, . . . ] denote the class of n-tuples equal to it up 
to cyclic permutations. We call such objects cyclic n- 
tuplcs, and its elements are naturally indexed by Z„. A 
closed walk of length n on a graph T is a cyclic n-tuple 
of oriented edges 

w = [e , . . . ,e„_i], e, € E, (80) 



such that It{ei) = / s (ej+i) for every i € Z„. The idea 
is that, given a graph, we can attach to it n-gons along 
closed walks. Note that the attachment can have two 
orientations, since given a closed walk w one could take 
the inverse walk w^ 1 := [e~_ 1; . . . , e^ 1 ] to describe the 
same attachment. Our definition of walks excludes the 
possibility of attaching the boundary of a disc along a 
walk consisting of a single vertex, something very useful 
in other contexts but not for our purposes. 

Let Wr denote the set of closed walks on the oriented 
graph r. An oriented 2-complex £ = (V, E, F, I s , I t , B) 
has the structure of a graph T = (V, E, I s ,It) plus a finite 
set F of faces and a boundary function B : F — > Wr- 
Just as we did for edges, we can consider the set F" 1 
of inverse faces setting B(f~ 1 ) = B(f)~ 1 . We also set 
F := FUF -1 . The discussion above explains how a topo- 
logical space M is related to this combinatorial structure 
S, and we will say that £ represents M and use them 
almost indistinguishably. In any case, our application 
to quantum error correcting codes only depends on the 
combinatorial point of view. Some examples will illus- 
trate the concept of 2-complex (see figure 11): 

• The sphere S. Take two vertices Vo,Vi, an edge 
e pointing from v to v\ and a face / with the 
boundary [e, e ]. 

• The projective plane P. Only a single vertex v, 
a single edge e and a single face / with boundary 
[e, e] are needed. 

• The torus T. This can be constructed with a ver- 
tex v, two edges e±, e2 and a face / with boundary 
[ei,e 2 ,e 1 " 1 ,e 2 " 1 ]- 

For any 2-complex X the Euler characteristic is 

XP):=\V\-\E\ + \F\. (81) 

S is said to be connected if its graph T is connected. 
£' = (V',E',F',r s ,I^B') is said to be a subcomplex 
of S if V C V, E' c E, F' c F, I' s c I s , I[ C I t 
and B' C B. As usual, we call components the maximal 
connected subcomplexes of E. Although we have defined 
X and connectedness in terms of S, they only depend 
upon the underlying topology. The same is true for Hi; 
its definition is our next goal. 

Consider a 2-complex E. For the shake of simplicity, let 
us introduce the notation A := V, Ai := E, A 2 := F. 
Let also the sets of 0-,l- and 2-chains be denoted Cj(£) 
with i = 0, 1, 2. They contain formal sums of elements of 
Ai with integer coefficients. We adopt the same conven- 
tions as for 0- and 1-chains for graphs. As in that case, 
Cj(£) ~ Z' Ai ' and A, is a natural basis of Cj(£). 

We introduce the boundary homomorphisms di : 
Ai — ► Aj_i for i = 1, 2. It is enough to give their value 
on a set of generators. We have: 

Vee£, 0i (e) = 7 t (e) - I s (e); (82) 
V/ef, d 2 (f) = c B{f y, (83) 
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FIG. 11: Planar representations of several 2-complexes. They 
have the topology of the sphere S, the projective plane P and 
the torus T. The identifiers for vertices and edges are the 
same as in the text. 




(a) (b) 

FIG. 12: A pair of 2-complexes embedded in R 3 . They rep- 
resent (a) the sphere S and (b) the torus T. 



where for any w = [el 1 , . . . , e° n ] € Wr, e% € E, Oi = ±1, 
we define 

n 

c w ■= ^2<Jzei, (84) 
i=\ 

Whenever the index i in di can be inferred from the con- 
text, we will omit it. A simple but fundamental property 
is 

d 2 = 0. (85) 

Let Zi(Yj) := kerdi, -Bi(E) := rancV The elements of 
Z\ are called cycles and the elements of B\ boundaries. 
We already encountered cycles in our study of the ho- 
mology of a graph. Note that B\ C Z\. Thus we can 
define 

Hi(S):=Zi(E)/Si(S). (86) 

Two cycles which represent the same element of the ho- 
mology group are said to be homologous. Boundaries are 
homologous to zero. If E consists of several components 
Ej, we have: 

ff 1 (E)~0^ 1 (EO. (87) 

i 

Our next goal is the definition of the first cohomology 
group H 1 ^). For i = 0, 1,2, let C J (E) denote the dual 
space of Cj(E), that is, 

C7*(E) :=hom(Ci(£),Z). (88) 



The elements of these spaces are called z-cochains. They 
can be regarded as the additive group of functions / : 
A; -> Z. Given a e A;, we define a* e C^E) by 

a*(a') = 6 aa ,, (89) 

where a' G A*. The set { ct* | cr e A 4 } is a basis of C J (E). 
For c* e C^E), c, e Cj(E), we let (cV;) := c l ( Ci ). 

For z = 1,2, we define the coboundary maps Si : 
C l_1 (E) — » C l (£) to be the dual homomorphism of <9j, 
that is, for every c 1 " 1 € C i_1 (E) and c l e Cj(E) we have 
((5c' _1 ,Cj) := (c J_1 ,9ci). Clearly, again omitting indices, 

S 2 = 0. (90) 

The set of cocycles and coboundaries are respectively 
Z X (E) := ker8 2 , B 1 ^) := ranJi. The first cohomology 
group is 

ff^E) := Z 1 (E)/S 1 (E). (91) 

Since they will be of interest when studying homolog- 
ical quantum error correcting codes, we collect here the 
following dual pair of properties. For any c\ <G C\ and 

yveV (Sv*,d)=0 <*=>■ CieZi; (92) 

V/£F (c\df) = c 1 ^ 1 ; (93) 

Also 

VueV,V/eF, (8v*,df) = 0. (94) 

We say that (<5v*, •) is a star operator and that (•, df) is 
a boundary operator, reflecting their geometrical nature. 
The name of the boundary operator is clear enough, but 
perhaps the star operator deserves some explanation. Let 
the star of a vertex v be the set 

star(u) := { (e, cr) e E x {1, -1} | I t {e a ) =v}. (95) 

Then we have 

sv* = ae *- ( 96 ) 

(e,(r)Gstar(v) 

D. Surfaces 

For a surface we understand a compact connected 
2-dimensional manifold. We already encountered sev- 
eral examples of surfaces constructed with 2-complexes, 
namely S, T and P. It is a fundamental result of sur- 
face topology that every other surface can be obtained by 
combination of these three; let us explain what is meant 
here by combination. 

Consider two surfaces, Mi and M 2 . Let Di, i = 1,2, 
be a subset of Mj homeomorphic to a closed disc and let 
its boundary be dDi. Let h : dDi — > dD 2 be a homeo- 
morphism. The connected sum of Mi and M 2 , denoted 
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FIG. 13: The connected sum of 2 tori to give a 2-torus in 
three steps: cutting, gluing and smoothing out. 




FIG. 14: The Moebius band. 
Mi jj M 2 , is defined as the quotient space of the disjoint 

o o 

union (Mi — D\) U (M 2 — D 2 ) under the identifications 

o 

x ~ h(x) for x £ dD\. Here Di denotes the interior of 
Di. Mi jjM 2 is a surface, and its homeomorphism class 
depends only upon the homeomorphism classes of Mi 
and M 2 . To gain intuition, figure 13 shows a connected 
sum of two tori to give a 2-torus. 

Let the Moebius band B be the topological space ob- 
tained as the quotient space of [0,1] x [0, 1] under the 
identifications [0, a;] ~ [1, 1 — x\. For a picture see figure 
14. A surface is said to be orientable if it does not con- 
tain a subset homeomorphic to B. A surface is embed- 
dable without self-intersections in R 3 iff it is orientable. 
S and T are orientable, but P is not. Define recursively 
gP := (g - l)P(jP for g > 1 and IP := P. Let also 
gT := (g - l)Tj}T for g > and 0T := S. No two of 
them are homeomorphic. gP is the sphere with n cross- 
caps and gT is the sphere with g handles or g-torus. gT 
and gP are said to have genus g. 

Proposition III.l Any orientable surface is homeomor- 
phic to gT for some integer n > 0. Any non- orientable 
surface is homeomorphic to gP for some integer n > 1. 

See for example [59] for a proof. We already presented 
above the standard 2-complexes representing P and T. 
gP can be represented by the 2-complex consisting of 
a vertex v, n edges {ai,...,a n } and a face / with 
B(f) = [ai, a\, . . . , a n , a n ]. gT can be constructed with a 



vertex v, 2n edges {ai, b\, . . . , a n , b n )} and a face / with 
B(f) = [ai,6i,af 1 ,&J" 1 ,...,a n ,6 n ,a- 1 ,6- 1 ]. Note that 
x(gT) = 2(1 — g) and x(j)P) = 2 — g. The correspond- 
ing homology and cohomology groups are Hi(gT) ~ 
H\gT) ~ T?9 and H^gP) ~ if^P) ~ Z^ 1 © Z 2 . 
The subgroup Z 2 appearing in the first homology group 
of non-orientable surfaces is called the torsion subgroup. 
It will play an important role when homological quan- 
tum error correcting codes for qudits of dimension greater 
than 2 are considered. 

Consider a topological graph T embedded in a surface 
M, that is, a homeomorphism between T and a subset of 
M. When M — T is a union of discs, we say that the em- 
bedding is a cell embedding. Clearly, such an embedding 
leads to a 2-complex whose faces are the mentioned discs. 
This raises the question of how to characterize combina- 
torially wether a 2-complex £ represents a surface or not. 
It is enough to give a condition such that for each vertex 
the corresponding point for the represented topological 
space has a neighborhood isomorphic to a disc. Let us 
first define the index of face / £ F on the 'corner' de- 
scribed by the ordered pair (e,e'), where e, e' £ E and 
It{e) = I s (e!). In plain words, the index counts the num- 
ber of times that the walk B(f) goes across the corner 
(ei, e 2 ). Formally, let B(f) = [e , . . . , e k -i] and 



|{ i £ Z fe | e = e; and e' = e i+1 }|; 



Then the index of / in (e, e') is 



indcx(/, e, e') 




ife- 1 



if e- 1 ^ e'; 



e . 



(97) 
(98) 



(99) 



So let v £ V be a vertex and let k := |star(v)|. We say 
that v is a surface vertex it there exists a cyclic fc-tuple 
S(v) = K°, . . . , e^- 1 ] such that star(v) = {{e u a,)}^ 1 
and 



if k= 1; 

^ 17 ' ' ~1l ifjfe>2,i-i = ±l; 

in other case, 

(100) 

where = is equality modulo k. Figure 15 illustrates the 
concept. Then, as a definition, a surface 2-complex is a 
connected 2-complcx such that all its vertices are surface 
vertices. We also need a way to distinguish orientability. 
We say that a surface 2-complex is oriented if 



0. 



(101) 



A surface 2-complex for which there is a suitable sign 
selection for faces so that it is oriented is said to be ori- 
entable. Figure 16 clarifies this definition. 

An interesting notion that emerges when considering 
the cell embedding of a graph in a surface is that of du- 
ality. The germ of this idea can be traced back to the 
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k= 1 k = 2 k=4 



FIG. 15: An illustration of the definition of surface ver- 
tex and the expressions (99) and (100). We have for ex- 
ample index(/i, ei, e^ 1 ) = 1, index(/2, ei, e?) = 1 and 
index(/4, 1 , e±) = 1. 




(a) (b) 



FIG. 16: Two planar representations of 2-complexes for (a) 
the sphere S and (b) the projective plane P. The identified 
vertices and edges are the same as in figure 11. S is shown 
with its faces oriented. On the other hand, P is not orientable. 
The picture shows an attempt to give a coherent orientation 
and the failure. Note that the impossibility comes from the 
presence of a Moebius band. It consists of the faces already 
oriented and the one with the interrogation sign. 

five regular platonic solids. Each of these polyhedra has 
a dual polyhedron whose vertices are the center points of 
the given one. For example, the tetrahedron is self-dual 
and the cube and the octahedron are dual of each other. 
The idea can be generalized. Given a cell embedding of 
a graph T in the surface M, the dual embedded graph T* 
is constructed as follows. For each face / a point /* is 
chosen to serve as a vertex for the new graph. For each 
edge e lying on the boundary of the faces /i and /2, the 
edge e* connects /* and crossing e once but no other 
edge or dual edge. Figure 17 shows a pair of examples. 

We now work out duality in the context of surface 
2-complexes. Consider an oriented surface 2-complcx 
£ = (V, E, F, I s ,I t ,B). We construct the dual 2-complcx 
S* = (V = F*,E' = E*,F' = V*,I' s ,I' t ,B') where 
V* := {v*\v £ V} and so on. There is a unique 
/ £ F such that (e*,df) = 1 (respectively -1) and 
we set I' s (e) = f* (respectively V t (e) = /*). For each 
v £ V, let S(v) = [e^°,...,el k _~i] be the cyclic k- 
tuple from the definition of surface 2-complexes. Then 
B'(v*) = [e^°,...,e* k _ 1 ' Tk - 1 ]. Now let the operator d 
take to «*, e to e*, and / to /*. Extend d linearly to 
act on any chain. Now, if we denote d* and S* the d and 
5 operators for S* , we have 

5*d = dd, dd* = Sd, (102) 

where the domains must be defined in the apparent way 




(a) (b) 

FIG. 17: (a) The standard cell embedding in the torus. Thin 
lines represent the dual graph. Notice self-duality, (b) A more 
complicated graph embedded in the torus and its dual. The 
torus is recovered from the plane model by identification of 
opposite edges. 

so that the composed function is well-defined. Finally, we 
observe that S* is oriented and X** ~ E. If one wants 
to extend the notion of duality to non-orientable surface 
2-complexes, Z 2 homology must be considered in order 
to eliminate orientation-related problems. We shall not 
dwell upon this here, however. 

Let us enlarge a bit the concept of surface. Take a 
surface M and a finite collection of disjoint sets {-Dj} 
such that each of them is homeomorphic to a disc. We 

o 

say that M — \J i Di is a surface with boundary. We al- 
ready encountered an important example of such an ob- 
ject, namely the Moebius band B. If one attaches to B 
a disc identifying homeomorphically its boundary with 
the rim of B, the projective plane P is obtained. We 
again need a combinatorial definition. Let S be a sur- 
face 2-complex and F' C F a collection of faces with no 
edge or vertex in common along the boundary walk. We 
say that S' := (V, E, F — F',I s ,I t ,B) is a surface with 
boundary 2-complcx. It is quite tempting to atempt an 
extension of duality to these broader class of 2-complexes. 
As the dual of a face is a vertex, it is apparent that the 
dual of a surface with boundary would be a 'surface with 
missing points'. Such an object is not a 2-complex, how- 
ever. To overcome this difficulty, relative homology can 
be considered. The relative homology of a complex re- 
spect to certain subcomplex is a topic in which we shall 
not enter, but it is worth mentioning that it would be 
perfectly suited to the error correcting code construction. 
Another possibility, is to construct the dual of a surface 
with boundary by identifying the corrcspondending ver- 
tices instead of deleting them. This construction leads us 
to what is called a pseudo-surface, a 'surface' which fails 
to be such a thing only in a finite set of points. From an 
homological point of view, the result is equivalent. See 
figure 18. 

For us the most important example of surface with 
boundary will be the ft.-holed disc Dh, h > 1. As a 
2-complex, Dh can be constructed with h + 1 vertices, 
2h + 1 edges and 1 face. Instead of giving explicitly the 
construction, we prefer to illustrate it with an example 
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(a) (b) (c) 

FIG. 18: (a) A 2-complex representing the disc with a hole, 
D\. (b) The dual of D\ in the form of a sphere with two 
vertices missing, (c) Another possibility for the dual of Di, 
obtained by identifying the missing vertices of the previous 
sphere. 




FIG. 19: A 2-complex representing the disc with four holes. 

in figure 19. We have H 1 (D h ) ~ H\D h ) ~ Z h and 
x(Dh) = 1 — h. The point of these perforated discs is 
that they have a nontrivial homology while still being a 
subset of the plane, something that we will find useful 
when physics come into play. 

E. Quantum homological codes 

From this point on we will be working with qudits of 
fixed dimension D. Unless otherwise stated, the homol- 
ogy considered will be always homology with coefficients 
in Z D . 

Before introducing homological quantum error correct- 
ing codes we still need a pair of definitions. Given a 2- 
complex E, let E — {ei}\Ji- Consider the isomorphisms 
hi : Ci(S) -» Z^ 1 and h 2 : C^E) -» Z 1 ^ 1 defined by 

\E\ 

hi(^A jei ) := (Ao,Ai,...,A| B |); (103) 

i=l 

\E\ 

h 2 (^A i e*):=(A ,A 1 ,...,A| B |). (104) 

i=l 

Let h : Ci(S) U C^E) -» Z^ E| be 

VcieCi(E) h(ci):=(0hi(ci)); (105) 

V^eC^E) h(c x ) := (ha^jO). (106) 

Then 

h^c 1 ) • hi(ci) = Mc 1 )* fih(ci) = (c\ci). (107) 



It is natural to use the notation u Cl := 0h( Cl ) and cr c i := 
(Th( c i) so that 

a c ia Cl = ( / 3((c 1 ,ci))c7 Cl c7 c i. (108) 

As we did for graphs, we can pull back the weight func- 
tion through hi and h 2 . Then we let the distance d(E) 
be the minimal weight among the representatives of non- 
trivial elements of Hi and H 1 . 

Theorem III. 2 Let E be a connected 2-complex. If 

F:=h[B 1 (E)]eh[B 1 (E)] (109) 

is generated by a Li. set, then setting Vc — V a symplec- 
tic [[n,k,d\] quantum error correcting code C is obtained 
with n = \E\, #!(£) ~ H 1 ^) ~ Z k D and d = d(E). 

Proof. The isotropy of Vc follows from (94) and (107). 
Also, from (92) we get hi[Zi(E)] = h 2 [B 1 (j:)}- L and 
h 2 [Z 1 {E)} =h 1 [ J B 1 (E)] ± , that is, 

V c =h[Z 1 (E)]®h[Z 1 (E)]. (110) 

But dimVc — dimVc = 2fc, and since dimhi[Zi(E)] + 
dimhalB^E)] = dim hi [Si (E)] + dim h 2 [Z 1 (E)] we get 
as desired H x (E) ~ H 1 (E) ~ Z k D . □ 

The condition that T be connected is just to avoid hav- 
ing a code which can be decomposed into two more simple 
ones. As for graphs, there is no point at all in considering 
disconnected 2-complexes; given such a disconnected 2- 
complcx with components Ej one can consider the wedge 
product of them, giving raise to the same code. 

The wedge product is obtained by choosing one vertex 
from each component and identifying them all. 

Because of the condition stating that the subspace 
(109) must be a linear subspace with a basis which is 
a linearly independent set in Z^ 1 , not every 2-complex 
can be used to produce codes for general qudits. For 
example, consider the case D = 4 in P, the projective 
plane. In this case Hi ~ Z 2 and thus a code cannot be 
constructed. The origin of the problem is in the torsion 
subgroup appearing in non-orientable surfaces. However, 
we can get rid of it if we only consider the case D = 2 
in this surfaces, as we shall do. Under this assumption 
and restricting attention to surface 2-complexes, we can 
give a more geometrical definition for the distance. Let S 
be a surface 2-complex, and let T be its graph. Let also 
Cy'(r) be the set of simple subcycles of T not homolo- 
gous to a point, and d'(E) the minimal length among the 
elements of Cy'(r). Then 

d(S) = minLi'(E), d'(E*)}. (Ill) 

To gain intuition on the construction of the codes, con- 
sider the special case of a graph as a 2-complex. In this 
case we obtain a pseudo-classical code, capable of cor- 
recting errors of the form <t x o whenever x is correctable 
in the corresponding classical code. 
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FIG. 20: A 2-complex giving rise to a [[25,1,5]] code. 




It is possible to construct homological quantum codes 
inspired by classical ones. Consider for example the 
graphs Cd, related to [d, l,d] classical linear codes. Join- 
ing d copies of Cd along vertices and attaching d(d — 1) 
faces, as shown in figure 20, gives a [[d 2 ,l,d]] code. In 
particular, for d = 3 we get Shor's original [[9, 1, 3]] code. 
Unfortunately, lim^oo — = 0, which is very different to 
the classical case. This fist example already shows that 
the length of quantum homological codes does not seem 
to behave very well when the distance grows. However, 
below we show that this is not the case when k grows. 

In general, homological quantum codes can be degen- 
erate. It is enough to have a vertex lying in less than d/2 
edges or a boundary with less than d/2 edges to have de- 
generacy. Such examples of degenerate codes will show 
up in the next section. 

F. Surface codes 

In this section we study homological quantum codes 
derived from 2-complexes representing surfaces. Such 
2-complexes are usually regarded as cell embeddings of 
graphs on surfaces, and so we will tend to use the lan- 
guage of topological graph theory. Note that the genus is 
directly related to the number of encoded qudits; codes 
derived from gT encode k = 2g qudits, and codes derived 
from gP encode k = g qubits. This can be put altogether 
using the Euler characteristic; cell embeddings of graphs 
on a surface M will give codes with 

k = 2- X (M). (112) 

As a first example of a surface code, figure 21 shows a 
self-dual embedding on P giving a [[9, 1, 3]] code. 

The whole problem of constructing good codes related 
to a certain surface relies on finding embeddings of graphs 
in such a way that both the embedded graph and its dual 
have a big distance whereas the number of edges keeps 
as small as possible. But let us be more accurate. 

Definition III. 3 Given a surface M and a positive inte- 
ger d we let the quantity /i(M, d) be the minimum number 
of edges among the embeddings of graphs in M giving a 
code of distance d. 

Since we do not know how to calculate the value of 
the function /x, we shall investigate some properties of 



FIG. 21: A self dual cell embedding in P. The projective 
plane is recovered by identifying opposite edges of the circum- 
ference. This embedding leads to a [[9, 1, 3]] code for qubits. 

this function. The problem of locality suggests also the 
introduction of a refinement of fx; the quantity fii(M,d) 
is defined as /i(M, d) but with the restriction that the 
graphs can have faces with at most I edges and vertices 
lying on at most / edges. Locality here means that we 
want that the vertex o$ v and face agf operators act on 
at most I qudits. Having operators as local as possible 
simplifies the error correction stage. We shall return on 
this issue below. 

We stress that in the case of non-orientable surfaces 
we only consider Z 2 homology. Keeping this in mind, we 
can state: 

Theorem III. 4 The function /x(M, d) is subadditive in 
its first argument, in the sense that given two surfaces 
Mi and M 2 

n{M4M 2 ,d) < f i(M 1 ,d)+n{M 2 ,d). (113) 

The proof is given in appendix B. 

The most simple orientable surface with nontrivial first 
homology group is the torus. In [19], a family of so called 
toric codes was presented, in the form of self-dual regular 
lattices on the torus. An investigation on other regular 
lattices on the torus led us to another system of lattices 
that demand half the number of qudits whereas it keeps 
the same good properties as the the first one; in particu- 
lar, vertex and face operators act on four qudits. In fact, 
in [19] only qubits were considered. Examples of both 
systems of lattices are depicted in figure 22, were the 
torus is represented as a quotient of the plane through 
a tessellation. In appendix C we show the optimality of 
our system. The original toric codes lead to a family of 
[[2d 2 ,2,d]] codes. Our lattices give [[d 2 + 1,2, d]] codes. 
This already shows that 

n(T,d) < in{T,d) < d 2 + 1. (114) 

Invoking subadditivity, we learn that fJ,(gT, d) is 0(x 2 ) in 
its second argumet, that is, it grows at most cuadratically 
with d. 

A closer examination of figure 22 reveals that the lat- 
tice giving a [10, 2, 3] code is a self-dual embedding of 
K 5 . This suggests considering self-dual embeddings of 
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FIG. 22: A graphical comparison between the original toric 
codes and their optimized versions. Thick lines are the border 
of tesserae and the torus is recovered as a quotient of the plane 
and its tessellation. (a),(b): The toric codes introduced in 
[19], for d — 3 and d = 5. (c),(d): The optimal regular toric 
codes for d — 3 and d = 5. 
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FIG. 23: The self-dual embedding of Kg in the 10-torus can 
be described using addition in Zg. It is enough to label the 
8-sided faces with i = 0, . . . , 8 and glue them altogether as 
the picture indicates. 




FIG. 24: The rate k/n vs. t/n for the codes derived from 
self-dual embeddings of the complete graphs K 4 i +1 (A) and 
for the optimized toric codes (□); O corresponds to the em- 
bedding of K$ in the torus. The quantum Hamming bound 
is displayed as a reference (solid line). 



K s , since such an embedding would give a [[(2)) (2) — 
2(s — 1),3]] code. In fact, these embeddings are possible 
in orientable surfaces with the suitable genus as long as 
s = l(mod 4) [60] and this family of codes with self-dual 
embeddings of complete graphs is enough to show that 
the coding rate k/n behaves as 



lim — 

k^oo n 



v 2 ff 
urn . — — 

s-oo v{gT, 3) 



1. 



(115) 



In order to verify this, note first that [i{gT, d) > 
fJ,(gT,l) = 2g. Let K(s) = (*) - 2(s - 1).' Due to 
subadditivity, for K(s)/2 < g < K(s + l)/2 we have 
p(gT, 3) < M^T, 3) + (g - ^)n(T, 3) = («) + 0(a). 

The limit (115) shows that the ratio k/n is asymptoti- 
cally one, and thus good codes can be constructed using 
surfaces. Figure 24 displays the rates for this family of 
codes and also for the optimized toric codes. The dif- 
ferences with figure 7 are apparent. However, the codes 
in the quantum case could be non-optimal, and thus the 
results are inconclusive. 



G. Planar codes 

We now focus on homological quantum codes derived 
from 2-complcxes representing surfaces with boundary. 
The situation is similar to the previous section and again 
we talk about cell embeddings of graphs. Note that for 
such a cell embedding of a graph on a surface with bound- 
ary, the boundaries are a subset of the graph. 

Surfaces with boundary offer more possible topologies 
to encode the same amount of qudits. If we remove from 
a g-torus I non-adjacent faces, H\ is enlarged with I — 
1 dimensions; removing a single face is useless since its 
boundary is a linear combination of the boundaries of 
the remaining faces. The non-orientable case is similar, 
because we only consider Z 2 homology. The results can 
again be collected using the Euler characteristic; given a 
surface with boundary M, not a surface, cell embeddings 
of graphs on it will give codes with 



k = l- x (M). 



(116) 



It is time to return on the issue of locality. Although 
topological codes are local, one has to face the problem 
of constructing a physical system with the shape of the 
surface on which the code lyes. At this point, the prob- 
lem of non-planarity arises; surfaces with non-trivial first 
homology group are not a subset of the plane, and so are 
difficult to realize experimentally. Among surfaces with 
boundary, however, there is such a planar family: the 
discs with h holes, D} Xl which encode h qudits. Figure 25 
displays the shape of non-correctable errors in Dh- 

An interesting point is that cell embeddings in gT giv- 
ing codes of distance d can be transformed to obtain cell 
embeddings in £>2 fl _i. The idea is to cut each of the han- 
dles of the torus, as shown if figure 26. The cut must be 
performed along a simple cycle of the graph, and so the 
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FIG. 25: Examples of non-correctable errors in D4. For clar- 
ity, the embedded graph is not shown. Thick lines represent 
typical elements of Hi , that is, cycles in the direct graph not 
homologous to zero. Dashed lines are elements of H 1 , in the 
form of cycles in the dual graph. 




FIG. 26: How to cut a torus of genus 2 to obtain a surface 
with boundary homeomorphic to a disc with 3 holes. 



edges of the cycle are duplicated in the process. These 
means that each cut introduces at least d new edges in 
the graph. On the other hand, the whole procedure pro- 
duces the lost of a single encoded qudit. A fundamental 
drawback of this method is that cocycles of length less 
than d could appear, thus diminishing the distance of 
the code. In such a case some additional edges could be 
added. However, it is also very possible that some edges 
become unnecessary after the cut: figure 27 shows an 
example. 

Another possible drawback of the cutting procedure is 
that the resulting embedding could be quite odd-shaped, 
and thus perhaps not very useful when true locality is 
necessary. In any case, one can always switch to more 
regular embeddings if the number of edges is unimpor- 
tant. Figure 28 displays such an embedding. 

It is possible to remove the condition that all faces 
must be homeomorphic to discs. In that case we are 
not dealing anymore with homology, but errors can still 
be visualized in a similar fashion. For example, Shor's 
[9,1,3] is displayed in figure 29. 




(a) (b) 

FIG. 27: (a) The result of cutting the self-dual embedding of 
K$ in T. (b) Some of the edges of the previous embedding 
can be deleted and still obtain a code of distance 3. 



FIG. 28: An example of an embedding in D4. The corre- 
sponding code has distance 3 and encodes 4 qudits. It is 
not difficult to generalize this embedding for general d and k; 
asymptotically the resulting code has n <x kd 2 . As usual, the 
growth is quadratic in d and linear in k. 




FIG. 29: A visualization of Shor's [9,1,3] code. It can also be 
considered as an embedding in the torus, since one face can 
be added for free. 



IV. CONCLUSIONS 

Quantum topology holds the promise of providing a 
mechanism for self-correcting errors without having to 
resort to constantly monitoring a quantum memory for 
error syndrome and error fixing. In this fashion, the func- 
tioning of a quantum memory would very much resemble 
the robustness of its classical counterpart. This is the 
main reason why it is very important to study quantum 
error correcting codes from a quantum topological point 
of view. In this paper we have acomplished this task 
by developing theorems characterizing homological quan- 
tum codes for qudits of arbitrary dimension D based on 
graphs embedded in surfaces of arbitrary topology, cither 
with or without boundaries, orientable or non-orientablc. 
Oricntability becomes an issue when trying to construct 
homological quantum codes using qudits of dimension 
D > 3 , due to the existence of a non-trivial torsion sub- 
group in the homology group. 

In doing so, we have realized that homological codes 
can also be well-defined in the classical case. This is in- 
teresting since not every classical code is of homological 
type. Nevertheless, we find that there exist a family of 
classical homological codes saturating the classical Ham- 
ming bound. 

As a result of our work, we have found that the prob- 
lem of constructing good quantum homological codes on 
arbitray surfaces relies on finding embeddings of graphs 
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in such a way that both the embedded graph and its dual 
graph have a big distance whereas the number of edges 
keeps as small as possible. This provides a connection be- 
tween the theory of quantum topological codes and topo- 
logical graph theory [60]. More specifically the problem 
of finding topological quantum codes is an instance of 
extremal graph theory which deals with the problem of 
finding maxima/minima of certain quantities defined on 
graphs. In our case, it is the distance of a quantum code 
wich has to be maximal on both the embedded graph 
and its dual. We have given an asymptotically optimal 
family of codes for the case of distance d = 3. We leave 
open the challenge of giving such optimal constructions 
for higher d. 
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APPENDIX A: GENERATORS OF Sp D (n) 

In order to proof that the homomorphism h introduced 
in (60) is onto, it is enough to exhibit a subset S C 
E Spain) such that h[S] generates Spoin). Consider 

• The Fourier operator on one qudit 

T:= vW)\k){l\\ (Al) 

k,iez D 

TXT ] = Z, TZT ] = X' 1 . (A2) 

• The operator on one qudit 

K:= £ f { l)U^+}iy m (A3) 

where the argument of <p must be evaluated in Z; 
KXK^ = f{l)XZ, KZK^ = Z. (A4) 

• The controlled NOT operator on two qudits 

tW:= Y, \k,l)(k,k + l\; (A5) 
kjez D 

UcvotX* ® Xiul Not =X l ® X i+ i , (A6) 

C/cNot^ ® ^'E^Not = zi ~ j ® 21 ■ (A7) 

The images under h of these operator on the first qudit (s) 
plus any qudit permutation generate Spoin) [56]. 



APPENDIX B: TOPOLOGICAL 
SUBADDITIVITY OF n 

We proof theorem III. 4. The assertion is quite trivial in 
the case d = 1. In order to proof it for d > 2, it is enough 
to construct an embedding of distance d in Mi$M 2 start- 
ing with two embeddings of distance d in Mi and M 2 in 
such a way that the number of edges does not increase; 
see figure 30. So let Si and E 2 be 2-complcxes of distance 
d representing respectively M\ and M 2 . We can suppose 
that neither of them is a sphere. Since d > 2, there exists 
an edge e\ in E^ which is not a self-loop. Let f\ be a 
face such that (/) = [uei, a, b, . . .], a e {1, —1}. We 
construct a new 2-complex T,[ introducing in Si a new 
edge e\ with the same source and target as e\ and chang- 
ing the boundary of f\ so that B^> (/1) = [ce'i, a,b, . . .]. 
We proceed in the same manner with £ 2 . Up to this 
point, we have performed the cutting step of figure 13 
and constructed two surfaces with boundary, T,[ and £ 2 . 
Then we construct E as a union of T,[ and £ 2 but identi- 
fying ei and e 2 in a single edge e, and similarly for their 
primed versions. Of course, the endpoints of e\ and e 2 
must be properly identified also, but the construction is 
clear enough so as to be self-explanatory. The resulting 2- 
complex is a surface, and that it represents the expected 
one follows from the two facts: it is orientable iff both 
Si and E 2 are orientable and x(E) = x(Ei) + x(E 2 ) — 2. 
We still have to check that its distance is d. The key 
observation is that e — e' is a boundary, in particular the 
boundary of the sum of all the faces in E' 1; properly ori- 
ented in the orientable case. Consider, for example, a 
simple cycle not homologous to zero that contains edges 
both from Ejv — {e, e'} and Ejv — {e, e'}; see figure 31. 
It must pass through each endpoint {v\, w 2 } of e exactly 
once. Then we can construct two simple cycles 71 and 7 2 
contained respectively in T[ and T' 2 . To this end we 'cut' 
7 in v\ and v 2 and glue again one of the pieces with e and 
the other with e'. At least one of the new simple cycles, 
say 71, is not homologous to zero in E, and thus in S^. 
Then its length is at least d, and the same is then true for 
the length of 7. Other possible simple cycles, including 
those in the dual graph, can be similarly worked out. 



APPENDIX C: OPTIMAL SELF-DUAL 
REGULAR TORIC CODES 

Let a cell embedding of a simplicial graph on a surface 
be a (v, /) regular cell embedding if the star of any vertex 
comprises v edges and the boundary of any face consists 
of / edges. On the torus, only the combinations (4,4), 
(3, 6) and (6, 3) are possible, since Euler's characteristic 
must be zero. We shall investigate here the self-dual case, 
(4, 4). In particular, given a distance d = 2i+ 1, we want 
to know which is the minimum number of edges in a (4, 4) 
regular cell-embedding on the torus such that its distance 
is d. 

We shall answer the question using homotopy. We 
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FIG. 30: The construction used to proof the subadditivity of 
(i. The first step is to perform a cut along a selected edge 
in each of the embeddings to be added. Then the resulting 
boundaries must be identified. 




FIG. 31: A simple cycle over the connected sum is divided in 
two, with each new simple cycle in one of the initial surfaces. 

say that an n-tuple w = (ei,...,e n ), G E, is a 
walk of length n if It(ei) — / s (e,+i), i = 1, . . . , n — 1. 
Its inverse is w^ 1 = (e^ 1 , . . . e^ 1 ). The empty walk is 
also a walk. If w = (. ..,e„) and «/ = (e' 1; ...) are 
such that It(e n ) = Zs(ei), then the composed walk is 



w + w' = (. . . , e„, e' 1; . . . ). If a walk is of the form 
w = wi + W2 + W3, and the boundary of a face (or its 
inverse) can be expressed as a walk as b = w 2 + W4, then 
we say that w and w' = Wi + w 4 1 + w 3 are homotopic 
and write w <~ w' . On a given embedding of a graph, 
we can choose any vertex v as a base point and consider 
the walks starting at v under the equivalence just stated. 
The resulting equivalence classes are the vertices of a new 
graph, naturally embedded in the universal cover of the 
surface under consideration. 

In the case of (4, 4) regular cell embeddings in the 
torus, the resulting graph is a infinite square lattice on 
the plane, as in figure 32. Let T be the original graph 
on the torus and T' the obtained graph on the plane. 



FIG. 32: An infinite square lattice on the plane. The vertices 
inside the dashed square are at most at distance 2 from the 
distinguished one. 



There is a natural projector p : V — > T taking vertices 
to vertices and edges to edges. Let v be the distinguished 
vertex in V representing the class of walks homotopic to 
a point. As in figure 32, we can consider the set of ver- 
tices at a distance at most t from v. If two of them have 
equal projections, say p{v\) — p(v 2 ), then there exists 
a walk going from v\ to v 2 of length less or equal to 2t 
such that its projection in T is not homotopic to a point. 
On a torus, this also means that it is not homologous to 
zero. Therefore, if T has distance d = 2t + 1, no such 
two vertices can exist. This means that F must have at 
least (d 2 + l)/2 vertices, and thus at least d 2 + 1 edges. 
As this minimal size is attained by the embeddings of 
section IIIF, we have the desired result. 
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